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Abstract 

Long's equation describes two dimensional stratified atmospheric flow over terrain which 
is represented by the geometry of the domain. The solutions of this equation over simple 
topography were investigated analytically and numerically by many authors. In this paper 
we derive a new terrain following formulation of this equation which incorporates the terrain 
as part of the differential equation rather than the geometry of the domain. This leads to 
new analytic insights about the solutions of this equation and enable us to compute steady 
state gravity wave patterns over complex topography. 



PACS 92.60.Gn, 92.60.Dj, 02.30.Ik 



1 Introduction 



Long's equation [1,2,3,4] models the flow of stratified incompressible fluid in two dimensions over 
terrain. When the base state of the flow (that is the unperturbed flow field far upstream) is 
without shear the numerical solutions (in the form of steady lee waves) of this equation over simple 
topography (i.e. one hill) were studied by many authors [5-13]. The most common approximation 
in these studies was to set Brunt- Vaisala frequency to a constant or a step function over the 
computational domain. Moreover the values of two physical parameters which appear in this 
equation were set to zero. (These parameters control the stratification and dispersive effects of the 
atmosphere - see Sec2.) In this (singular) limit the nonlinear terms and one of the leading second 
order derivatives in the equation drop out and the equation reduces to that of a linear harmonic 
oscillator over two dimensional domain. Careful studies [8] showed that these approximations set 
strong limitations on the validity of the derived solutions [9] . 

Long's equation also provides the theoretical framework for the analysis of experimental data 
[14,15,29] under the assumption of shearless base flow. (An assumption which, in general, is not 
supported by the data). An extensive hst of references appears in [16,17,18]. 

An analytic approach to the study of the solutions of this nonlinear equation was initiated 
recently by the current author [19,20,21]. We showed that for a base fiow without shear and under 
rather mild restrictions the nonlinear terms in the equation can be simplified. We also identified 
the "slow variable" that controls the nonlinear oscillations in this equation. Using phase averaging 
approximation we derived for self similar solutions of this equation a formula for the attenuation 
of the stream function perturbation with height. This result is generically related to the presence 
of the nonlinear terms in Long's equation. The impact that shear has on the generation and 
amplitude of gravity waves was investigated by us in [20]. A new representation of this equation 
in terms of the atmospheric density was derived in [21]. 

One of the weak aspects of Long's equation is related to the fact that the terrain is represented 
by the shape of the domain and the boundary conditions. As a result the impact of different terrains 
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on the solution of this equation can only be studied numerically. Furthermore discretization errors 
which occur in the representation of the terrain render it impractical to consider complex terrain. 
In part these errors are due to the scale of the terrain relative to the computational domain. 
Accordingly only simple topographies which were represented by one hill were considered in the 
literature. Furthermore even for these simple topographies only approximate boundary conditions 
were applied at the terrain. (See discussion in Sec 2). 

With this motivation it is our objective in this paper to derive a terrain following formulation 
of Long's equation in which the terrain is incorporated as part of the coefficients of the differential 
equation, and the computational domain is always a rectangle. This new representation makes it 
possible to derive new analytic insights about the solution of this equation in some limiting cases. 
It will make it easier also to study how the solution varies as a function of the terrain and other 
parameters that appear in the equation. 

The plan of the paper is as follows: Sec. 2 presents a short review of Long"s equation and some 
aspects of its solutions. In Sec. 3 we derive the new formulation of this equation. Sec 4 considers 
some analytic and geophysical aspects of this new formulation while Sec 5 compares its numerical 
solution over three different terrains. These simulations are motivated by recent experiments in 
the Alps region to educe properties of gravity waves from experimental data [14]. We end up in 
Sec 6. with a summary and conclusions. 

2 Long's Equation - A Short Overview 

In two dimensions (x, z) the flow of a steady inviscid and incompressible stratified fluid is modeled 
by the following equations: 




+ wpz = 



(2.2) 




(2.3) 
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p{uW:r + WWz) = -p^ - PQ (2.4) 

where subscripts indicate differentiation with respect to the indicated variable, u = {u, w) is the 
ffuid velocity, p is its density p is the pressure and g is the acceleration of gravity. 

We can non-dimensionalize these equations by introducing 



X _ A^o - u _ LNq 

X = Z = T^Z, u = — , W = ^^r, w 



L' '-JT,' "-f/o' ---Ul 
p _ No 



p = p = —rr—p (2.5) 

Po gUoPo 

where L represents a characteristic horizontal length, and Uq,Pq represent respectively the free 
stream velocity and averaged base density (i.e. here po is a constant). Nq is an averaged value of 
the Brunt- Vaisala frequency 

N' = -^^ (2.6) 

Po dz 

where po = Poi^) is the base density. 

In these new variables eqs fl2.ip - 02.41) take the following form (for brevity we drop the bars) 

u^ + w, = (2.7) 

up^ + wp^ = (2.8) 
(3p{uu^ + wu^) = -p^ (2.9) 
I3p{uw^ + ww^) = -p^'^iPz + p) (2.10) 

where 

13= Ml (2.11) 

9 

(3 is the Boussinesq parameter [13] (this name has nothing to do with the "Boussinesq approxi- 
mation") which controls stratification effects (assuming Uq ^ 0) and p is the long wave parameter 
which controls dispersive effects (or the deviation from the hydrostatic approximation). In the 
limit p = the hydrostatic approximation is fully satisfied, [10,11]. 
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In view of eq. (12.71) we can introduce a stream function \1/ so that 

u = ^„ w = -^^ . 



(2.13) 



Using this stream function we can rewrite eq. (12.81) as 

J{p,vI/} = (2.14) 

where for any two (smooth) functions /, g 

jlfg\ = ^^^^^ (215) 
dx dz dz dx 

Eq. (I2.14P imphes that the functions p, \l/ are dependent on each other and we can express each 
of them in terms of the other. Thus we can write \1/ as ^I/(p) (or p as p(^I/) [21]). 

After a long algebra one can derive the following equation for \E' [22,1,13] 



where 



z+^{^i+p'^i: 



S'(^) (2.16) 



N^m = (2.17) 

fjp 

is the nondimensional Brunt- Vaisala frequency. We observe that in this definition A^^ is a function 
of \E'. (As a result it can be an additional source of nonlinearity in eq. (12.161) ) This is in contrast 
to the previous definition of this quantity in eq. (12. 6p which depends only on the base state. In 
the following we assume without loss of generality that the direction of base flow is from left to 
right along the x-axis. Furthermore we assume it to be a function of z only. 

5'(\I') is some unknown function which is determined from the base flow. To carry out this 
determination of S we consider eq. (I2.16P as x — oo and express the left hand side of this 
equation in terms of \I' only. (Assuming that disturbances do not propagate far upstream [16,18]). 
Eq. (12.161) is referred to as Long's equation. 

For example if we let 

lim <if{x,z) = z (2.18) 



X— »— oo 
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i.e consider a shearless base flow with lim u{x, z) = 1 then 

X— >— oo 

Si^qj) = -N\^) (^^ + (2.19) 

and eq. (12.161) becomes: 

^zz + ^'^xx - iV'(^)[^ - ^ + ^(*^ + fi^^l - 1)] = 0. (2.20) 



It is evident from this derivation that different profiles for the base flow as x — — oo will lead to 
different forms of S{'^)[20]. 

For a general base flow in an unbounded domain over topography with shape f{x) and maxi- 
mum height Ho the following boundary conditions are imposed on \E' 

lim ^(a:,z) = ^o(^) (2.21) 

I— > — oo 

H N 

'^{x,Tf{x)) = constant, r = (2.22) 

Uq 

where the constant in eq. (12.22p is (usually) set to zero. As to the boundary condition at x — > oo 
it is appropriate to set 

lim \l'(x, z) = "^oiz) 

(in spite of the fact that Long's equation contains no dissipation terms). However over finite 
computational domain only radiation boundary conditions can be imposed in this limit. Similarly 
as z ^ oo it is customary to impose (following [7]) radiation boundary conditions. (The imposition 
of these boundary conditions is discussed in detail in S'ec 4.1). 

For the perturbation from the shearless base flow 

= ^ - 2 (2.23) 

eq. fl2.20p becomes 

r^.. - a'r]l + fi\vxx - a'vl) - N\r]){Prj, - v) = (2.24) 
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where 



Ar2 



(2.25) 
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We observe that when |r| ^ 1 the boundary condition fl2.22p can be approximated by 



r]{x,0) 



rf{x). 



(2.26) 



When is constant eq. (12.24p is invariant with respect to translations in x, z and hence admits 
self-similar solutions of the form r] = f{kx + mz) [19]. These solutions are interpreted as gravity 
waves that are generated by the flow over the topography. 

From a numerical point of view it is a common practice [7,8,13] to solve eq. (12.241) in the limit 
/3 = and /i = with constant over the domain. However observe that the definition of in 
Long equation is given by (I2.17P and it depends on In some other numerical simulations the 
computational domain is divided into subdomains where is constant in each subdomain but 
this led to numerical instabilities at the interface between these subdomains. 

In these limits Eq. (I2.24p reduces then to a linear equation 



We observe that the limit (3 = can be obtained either by letting Uq ^ or Nq ^ 0. In the 
following we assume that this limit is obtained a.s Uq (so that stratification persists in this 
limit and the leading term in A^q is not zero). 

Eq. (12.270 is a singular limit of Long's equation as one of the leading second order derivatives 
drops when = and the nonlinear terms drops out when (3 = and A^ is constant. This 
approximation and its limitations were considered numerically and analytically [6,7,19,20] and 
was found to be justified only under strong restrictions even under the assumption that the base 
flow is shearless. Nevertheless it is used routinely in the actual analysis of atmospheric data 



(2.27) 



[14,15,16]. 



The general solution of eq. (12.270 is 



Tji^x, z) = q{x) coii{Nz) + p{x) sin(A^2;) 



(2.28) 
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where the functions p{x), q{x) have to be determined so that the the boundary conditions derived 
from eq. fl2.22p . fl2.26p and the radiation boundary conditions are satisfied. These lead in general 
to an integral equation for p{x) and q{x) and it easy to show [13] that p{x) = H[q{x)] where 
H[q{x)] is the Hilbert transform of q{x). The boundary condition on the terrain becomes; 

q{x) cos{TNf{x)) + H[q{x)] sm{TNf{x)) = -Tf{x) . (2.29) 

This integral equation has to be solved numerically [6,7,13]. 



3 Terrain Following Formulation. 



To derive a terrain following formulation of Long's equation which incorporates the terrain in the 
coefficients of the differential equation (rather than the shape of the domain) we introduce Gal- 
Chen transformation. If the height of the (bottom) terrain is described by a sufficiently smooth 
function z = h{x) and the height of the computational flow region is flnite, i.e. h{x) < z < H, 
where H is a. constant, then this transformation is given by 

- h{x) 



X 



H 



H-h{x)' 



Under this transformation we have 

d_ 

dx 

where 



dx dz ' 



d_ 

dz 



H 



G 



12 



'G H-h{xy " 
Furthermore the expression of the Laplace operator becomes 



]_d_ 



- 1 h'ix) 



dx' 



+ 



G 



+ {G 



12\2 



d'^ d'^ 
dz?' dxdz 



+ 



dG'^ ^^,,dG'^ 



dx 



dz 



d_ 

dz 



(3.1) 



(3.2) 



(3.3) 



(3.4) 



Under this transformation the continuity equation (12.71) becomes 

1 dw 



du ^ _^ 
dx dz 



G dz 



0. 



(3.5) 
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However if we introduce 

v = ^{w + VGG^\) (3.6) 



then it is a simple algebra to show that eq. fl3.5p can be rewritten as 

+ 1(75.0 = 0. (3.7) 

From this equation we see that we can introduce a "terrain following stream function" ijj so that 

u = VGu = ^, v = ^v = -^. (3.8) 
oz ox 

Multiplying eq. (12. 8p by ^/G we can rewrite this equation in the following form: 
Using eq. (13.81) this can be rewritten as 

J{p,ij} = (3.10) 

where J is defined as in eq. (I2.15P but with differentiations with respect to {x, z). Equation (I3.10p 
implies that p{x, z) = p{i/^{x, z)) (and vice versa). 



To eliminate the pressure term from eqs. (12. 9p and (I2.10p we differentiate (12. 9p by z and apply 
the operator /^^^ to (I2.10p and subtract. We obtain 

I3p^p^{uw^ + ww^) - Ppz{uu^ + wuz) + (3p^p{uw^ + ww^)^ - (3p{uu^ + wu^)^ = -p^. (3.11) 

Using eq. (12.80 the first two terms in this equation can be written as 

(3p^p^{uw^ + ww^) - Ppz{uu^ + wu^) = (3 [p'^{-p^ww^ + p^ww^) - p^uu^ + p^.MM^] (3.12) 

/^r^2|2 2\ /2|2 2\1 \ f 2 , 2 2\ / 2 , 2 2\ 1 

= 77 [Px{u + p W )z- Pz[u + p W )^\ = j= [ps{u + p W )2- Pz[u + p W )s\ = 

^ 2v G 

^ Ir 2 I 2 2t 

J{p, u + p w \. 



2^/G 
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Using eqs (13. 6p and (13. 8p to re-express + fs^w"^ we have 



2\/G 2v Ct 



(3.13) 



The third and the fourth terms in eq. (13.111) can be rewritten using (12.71) as 

Pi? p{uw^+wWz)^-Pp{uu^ + wUz)z = Pp [u{p'^w^ - Uz) + v{p^w^ - Uz)] = — ^J{i/j,x} (3.14) 



G 



where x = l?Wx — Uz is the vorticity. Expressing x terms of il) we have 



(3.15) 



where 



dx 



dz 



dx"^ ' dxdz 
is the "terrain following Laplace operator". 

Finally for the right hand side of eq. (13.111) we have 



d_ 

dz 



- Px 



J {p. 9] 



where 



g{x, z) = z + h{x) ( 1 — 



G 



H 
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(3.17) 



Combining all the results contained in eqs. (]3.12p - (]3.17p we can re-express eq. (13. lip in the 
following form: 



p\i;,r + 2p'G''rx^, + ( ^ + p?{G''f ] {^,f 



N\i;)g{x,z) 



= 
(3.18) 



where N'^{il)) is defined as in eq. (I2.17p . Hence it follows that, 



p\^l^,f + 2i?G''i,,^-z + ( ^ + p\G''f ) {i,,f 



-N\i,)g{xrz) = S{',p). (3.19) 
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This is the terrain following form of Long's equation. At this juncture it might be asked why 
one can not "save" this derivation and apply the terrain following transformation (13. ip directly 
to f l2.16p . Doing so will yield an extremely complicated equation. This has been avoided in our 
derivation by the use of the "terrain following stream function" in (13. 8p . 



To determine the function S{iIj) in eq. (I3.19P we assume that 

lim h{x) = 

X— >— oo 

and that (as an example) ip satisfies 

lim ip{x, z) 

It follows then that 



z. 



and Long's equation becomes: 



(3.20) 



(3.21) 



g{x,z) - ^"2 



0. 



In this representation the flow domain is a rectangle [a, &] x [0, H] or an infinite stripe [— oo, oo] x 
[0, H]. The boundary condition at the bottom topography is 

u ■ n = 

where n is the normal to the topography which is described by the curve h{x). Hence this normal 
is given by n = {—h'{x), 1). Using (I3.6p . (l3.8p this leads to the boundary condition 



i{j{x, 0) = constant 



(3.22) 



and this constant can be chosen to be zero. The other boundary condition that has to be imposed 
on t/) is a radiation boundary condition at z = if (which implies that the outgoing wave is not 
reflected by the boundary). 
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To obtain an equation for the perturbation from the base state we set 



z + r]{x, z). 



(3.23) 



Substituting this in eq. fl3.2ip we obtain the following (exact) equation for rj 




4 Analytic solutions of Long's equation 

In the traditional representation of Long's equation the topography determines the shape of the 
flow domain and as a result it is not feasible to obtain analytic solutions to this equation even in 
some limits of the parameters fi and /i . We now show that this problem can be overcome in some 
limiting cases when the terrain following formulation of this equation is used. 

We consider two limiting cases /3 = 0, /i = and /3 7^ 0, = we also assume N'^{ip) = constant. 
For brevity we drop in the following the bars over x, z. 

4.1 The Limiting Case = 0, = 

In this case eq. f l3.2ip simplifies to 





whose general solution is 



ijj = A{x) cos(z/z) + B{x)sin{h'z) + z + h{x){l 



z 




Here v = N\fG and A{x), B{x) are functions which have to be determined from the boundary 
conditions . 
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The boundary condition fl3.22p implies A{x) = —h{x). To determine B{x) we must apply the 
radiation boundary condition as 2; — ^ oo on the solution. To this end we must insure that the 
vertical group velocity of the wave is positive. Using the dispersion relation for hydrostatic flow 
given in [16,p.l81] this group velocity is: 

c, = ^Ih^ (4.3) 

where k is the horizontal wave number. We deduce then that the vertical group velocity is positive 
when ku > 0. 

To impose this condition on the solution (14.21) we express A{x), B{x) in Fourier integral form 

/oo /"OO 
a{k)e'^''dk, B{x) = / h{k)e'^''dk (4.4) 
-00 J —00 

where k is the horizontal wave number. We deduce then that the solution (14.21) can be written as 

iP = ^\ [ {a{k) - ib{k))e'^'''^+''''>dk + [ {a{k) + ib{k))e'^'''^-'''Uk\ + \z + h{x){l - ^) 
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(4.5) 

To satisfy the radiation boundary condition for 2; —> 00 the first and second integral must vanish 
for A; < and A; > respectively. Therefore a{k) and h{k) must satisfy 

a{k) = -isgn{k)b{k) (4.6) 

which implies that B{x) is the Hilbert transform of A{x) = —h{x). This represents a complete 
analytic solution of Long's equation for this limiting case. 

We compare now these analytic results with the solution methodology that has been used 
previously in the literature as was discussed Sec 2. First we note that this analytic solution 
requires only the direct (and simple) computation of the Hilbert transform of the terrain function 
h{x) . This is a straightforward procedure even if it has to be done numerically. On the other 
hand to compute q{x) using (I2.29P requires in general the solution of an integral equation. To do 
so one must use an iterative algorithm which might turn out to be unstable or non-convergent 
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over complex terrain. Furthermore there is the issue of applying the boundary conditions on ifj at 
the terrain. To this end the procedure discussed in Sec 2 requires the use of the approximations 
that lead to (12.261) . As a result the equation that is used to compute q{x) (eq. f l2.29p ) is also an 
approximate equation which will yield at best approximate solution for this function. On the other 
hand the application of the boundary conditions using the procedure discussed in this section is 
exact and does not place constraints on the height of the terrain. 

From an experimental geophysical point of view it has been a common practice to assume that 
the gravity wave generated by a flow over terrain is of the form sai{kx + mz) (or similar) [14,15,28]. 
This has led to difficulties in the eduction of this wave from experimental data. Our results 
show that this form of the wave is incorrect (at least in principle). Furthermore the numerical 
simulations that we carry in the next section demonstrate that complex terrain can alter drastically 
the shape and amplitude of this wave due to interference effects. 



4.2 The Limiting Case ^ 0, = 



In this limiting case eq. (13.211) becomes 



^ 2 



+ z + h{x){i--: 



(4.7) 



Since this is a nonlinear equation we can find an approximate analytical solution using first 



order perturbation expansion under the assumption that 
practical situations). Expressing i/j approximately as 



<^ 1 (which is satisfied in most 



and substituting this expression in (14. 7p we obtain to order zero and one in the parameter the 
following equations 



dz'^ 



+ GN'iJo = GN' 



z + h{x){l 



H' 



dz^ 



l + GN^ij,-aH^? = -Ga^ 



dz 



(4.8) 
(4.9) 
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The boundary conditions on "00, "^i cire given by eq. (13.221) at z = and radiation boundary 
conditions as z ^ oo. 

Solving these (hnear) equations for ipQ and i/^i we obtain the following expressions for their 
solutions 



xjjQ = A{x) cos^uz) + B{x) sm(uz) + [z + h{x){l — — )] 

H 

ijji = Ci{x) cos{uz) + C2{x) sin(i^z) + z) + /2(x, z) + fsix) - ^ 



(4.10) 
(4.11) 



where 



fi{x,z) 



{A^{x) - B^{x))cos{2vz) + 2A{x)B{x)sin{2vz) 

6 

[H — h{x)){{uA{x)z + B{x)) cos(z/2;) + i/B{x)z sin^uz) 



_ A{xf + B{xf 1 



h{x) 



The determination of the functions A{x), B{x), Ci{x), C2{x) from the boundary conditions can be 
done using the same procedure outlined in the previous subsection. However as it is algebraically 
cumbersome we omit the details. 



5 Numerical Simulations over Complex Terrain. 

In the previous section we presented analytical solutions to Long's equation in some limiting 
cases. In general one has to resort to numerical simulations of Long's equation. However as it 
was said in the introduction the terrain following form of Long's equation obviate the need to 
discretize the terrain in these numerical procedures and therefore lead to better representations of 
complex terrains. Also the application of the boundary condition at the terrain is being simplified 
considerably and is treated exactly contrary to the approach discussed in Sec 2 (see eqs. (12.261) 
and fl2:29|) ). 

These simulations are carried in order to validate our equations and explore the different 
flow patterns (and possible special effects) that can be predicted over complex topography. The 
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geophysical motivation for these simulations is related to recent experiments to measure and educe 
gravity waves and their properties over the Alps from balloon data. This endeavor faced several 
difficulties. In part these difficulties can be traced to the fact that the terrain over which the 
measurements were made contained (at least) two summits rather than one [15]. 

In this section we carry simulations over three terrains. These consist of one, two and three 
hills with the following shape functions 

/(^) = (IT^' ^^-'^ 

^^^^ " (1 + x2)3/2 + (1 + (2: - 5)2)3/2 ' ^^-^^ 
^^^^ " (1 + x2)3/2 + (1 + - 5)2)3/2 + (1 + (a; - 10)2)3/2 " ^^'^^ 

In all cases we solved fl3.24p for the perturbation with the following parameters; 

r = 0.25, A^ = l, /5=1.10-^ ^ = 0.1. (5.4) 

Here represents the nondimensional Brunt- Vaisala frequency which was defined in ( 12.171) . 

To solve for the perturbation rj over a finite two dimensional domain which after the trans- 
formation (13.11) is represented by [a, c] x [0, H] we imposed at the boundary z = the condition 
r] = 0. Radiation boundary conditions were imposed at x = c and z = H. These are necessary to 
avoid refiection of the outgoing wave. To implement these boundary conditions we used "sponge 
boundaries" at x = c and z = H [23,24]. The sponge damping constant was calibrated to suppress 
wave refiection at these boundaries. In addition open boundary conditions were used at these 
boundaries themselves viz. we let 

^ic,z) = 0, ^^{x,H) = 0. 

In the simulations described below we let (in nondimensional units) a = —25, c = 25 and 
H = 25 with a grid of 801 x 401 points. The sponge layer (at a, c) consists of 20 grid points. 
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To solve (13.241) under these settings we used Matlab[28]. Central finite differences approxi- 
mations were used to discretize Long's equation on the domain grid and a fixed point iterative 
algorithm was implemented to solve the resulting equations. The convergence criteria for the iter- 
ations was that the step error \r]m+i — "^ml was less than 1.10~^ where m is the iteration number. 

Convergence was achieved in less than 100 iterations. 

— It should be noted however that 77 in these simulations represents the perturbation from the 
"terrain following stream function" which was defined in (13. 8p . Using (13. 6p . (13. 8p we find that the 
perturbation to the base flow field Up, Wp is given by 

, , h(x) H dr] . . . ^, 

Wp{x, z) = -||(x, -z) - G''{x, -z) + |^(x, -z)^ (5.6) 

where z is defined in (13. ip . From this vector field we can compute the "regular" stream function 
0(x, z) which is defined as 

(j){x, z) = z + vi^x^ z) (5.7) 

where z/(a;, z) satisfies 

= 91' = -Yx ^'-^^ 

A zoom-plot of (p{x, z) (on part of the computational domain) on the region around the three 
topographies discussed above is presented in Figsl,2,?). We observe that the phase lines of the 
stream function in all these figures tilts with height. In all three figures the gravity waves are in 
the lee of the topography. 

These figures demonstrate the impact that complex terrain can have on the structure of gravity 
waves. 

Fig 1 displays the results of the simulation for one hill (described by (15.10 ). The plot corre- 
sponds to the classical results that appeared in the literature for this case. The results for two hills 
in Fig 2 show some wave activity between the two hills and less organized waves with somewhat 
smaller amplitude in the lee of the two hills. This might be due to interference effects between 
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the waves generated by the two hills. For three hills the simulation shows the same wave activity 
between the two hills (on the left) as in Fig 2 but then a "quiet" zone between the second and 
third hills in which is almost constant. The waves in the lee of the three hill appear regular as 
in Fig 1. 



6 Summary and Conclusions. 

We derived in this paper a terrain following formulation of Long's equation in which the topography 
is "absorbed" in the coefficients of the differential equation representing the flow rather than being 
part of the boundary conditions. We used this representation to solve Long's equation analytically 
in some limiting cases and numerically over complex topography. The new formulation also opens 
the possibility to develop analytical estimates which compare the solutions of this equation over 
difi^erent topographies . 

Prom a geophysical point of view it well known that some present models for the generation 
of gravity waves over estimate this effect [25,27]. Partially, this is due to the fact that they use 
oversimplified representation of the terrain. Purthermore they do not take into account the effects 
that are due to complex terrain (as demonstrated by our simulations). We believe that the new 
form of Long's equation will make it easier to consider more realistic representations of the terrain 
and its effect on the generation and propagation of gravity waves. 
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List of Captions 



Fig. 1 Zoom-Contour plot of the stream function (p over one hill 
centered at x = with N = 1, /3 = 0.01, fi = 0.3, r = 0.25. 

Fig. 2 Same as Fig 1 with two hills centered at a; = and x = 5. 

Fig. 3 Same as Fig 1 with three hills centered at a; = —5, and 5. 
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Contour plot of the Stream Function 
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Contour plot of the Stream Function 
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